Draft version October 13, 2008 

Preprint typeset using I^T^X style cmulatcapj v. 11/26/04 



oo 
O 
O 
(N 
•*-> 
O 

O 

m 



43 

6 



> 
O 

cn 
o 
oo 
o 



X 



DETECTION OF X-RAY EMISSION FROM THE VERY OLD PULSAR J0108-1431 
G. G. Pavlov 1 , O. Kargaltsev 1 ' 2 , J. A. Wong 1 , and G. P. Garmire 1 

1 Dept. of Astronomy and Astrophysics, Pennsylvania State University, PA 16802; pavlov@astro.psu.edu 
2 Dept. of Astronomy, University of Florida, Bryant Space Science Center, Gainesville, FL 32611; oyklOO@astro.ufl.edu 

Draft version October 13, 2008 

ABSTRACT 

PSR J0108-1431 is a nearby, 170 Myr old, very faint radio pulsar near the "pulsar death line" in 
the P-P diagram. We observed the pulsar field with the Chandra X-ray Observatory and detected 



a point source (53 counts in a 30 ks exposure; energy flux (9 ± 2) x 10 5 ergs cm 2 s 



the 



0.3-8 keV band) close to the radio pulsar position. Based on the large X-ray/optical flux ratio at the 
X-ray source position, we conclude that the source is the X-ray counterpart of PSR J0108-1431. The 
pulsar spectrum can be described by a power-law model with photon index r w 2.2 and luminosity 
130 er S s s \ or by a blackbody model with the temperature kT w 0.28 keV 
and bolometric luminosity Lboi ~ 1-3 x 10 28 df 30 ergs s -1 , for a plausible hydrogen column density 
Yh = 7.3 x 10 19 cm~ 2 (di3o = rf/130pc). The pulsar converts ~ 0.4% of its spin-down power into 
the X-ray luminosity, i.e., its X-ray efficiency is higher than for most younger pulsars. From the 
comparison of the X-ray position with the previously measured radio positions, we estimated the 
pulsar proper motion of 0.2 arcsec yr _1 (V± ~ 130di3o km s _1 ), in the south-southeast direction. 

Subject headings: pulsars: individual (PSR J0108-1431) — stars: neutron — X-rays: stars 



1. INTRODUCTION 

Because of the observational selection, most of radio 
pulsars detected in X-rays are much younger and more 
powerful than "typical" ones, and they cannot be consid- 
ered as a representative sample of the whole radio pulsar 
population. Therefore, to understand the evolution of 
X-ray pulsar properties, it is important to expand the 
current sample toward lower spin-down powers, E, and 
larger spin-down ages, r. X-ray emission from such low- 
powered, old pulsars is expected to be composed of non- 
thermal (magnetospheric) emission and thermal emission 
from the heated polar caps around the neutron star (NS) 
magnetic poles. (Old NSs are too cold to observe thermal 
emission from the bulk of the NS surface in X-rays; see, 
e.g., Yakovlev & Pethick 2004.) The magnetospheric X- 
ray emission is routinely seen in both young and old pul- 
sars. It is characterized by a power-law (PL) spectrum, 
with a typical slope (photon index) 1 < T < 2, and rel- 
atively sharp pulsations. The X-ray emission from polar 
caps, heated by relativistic particles accelerated in the 
pulsar magnetosphere, is expected to be most efficient 
in old pulsars (e.g., Harding & Muslimov 2001, 2002). 
It should show a thermal spectrum, with a temperature 
kT ~ 0.1-0.3 keV, and smoother pulsations. Since both 
the magnetospheric and thermal components are pow- 
ered by the NS spin-down, their luminosities are expected 
to be proportional to the spin-down power, Lx = VxE, 
where rjx < 1 is the pulsar X-ray efficiency. 

As of February 2008, Chandra and XMM-Newton ob- 
servations of eight old, low-powered pulsars (lMyr < 
t < 20 Myr, 31.9 < logE < 33.6) have been reported, 
of which seven were detected (Becker et al. 2004, 2005, 
2006; Zavlin & Pavlov 2004; Zhang et al. 2005; Tepede- 
lenlioglu & Ogelman 2005; Kargaltsev et al. 2006; Hui 
& Becker 2007; Misanovic et al. 2008) 1 . Most of them 

1 After submitting this paper, Gil et al. (2008) reported XMM- 



showed rather faint X-ray emission, such that it was hard 
to accurately characterize their spectra and detect pul- 
sations. The spectra of all the detected pulsars can be 
satisfactorily fitted with a PL model, but the PL slopes, 
T ~ 2-4, are, on average, steeper than those found in 
young pulsars, and the hydrogen column densities, Yjj, 
often exceed those estimated from the pulsar dispersion 
measure at the usual assumption of 10% ionization of 
the interstellar medium (ISM). (An extreme example is 
PSR B1929+10, for which the Y H found from the PL fit 
corresponds to 0.4% of ISM ionization; Misanovic et al. 
2008). Some of the spectra can be fitted with a black- 
body (BB) model, but the emitting areas are usually 
much smaller than the conventional areas of polar caps. 
The spectra can be also described by a two-component, 
PL+BB, model, but usually with poorly constrained pa- 
rameters. The measured X-ray efficiencies of old pul- 
sars show a large scatter, rj ~ 10~ 4 -10~ 2 in the 1-10 
keV band, but, on average, they seem to be higher than 
those for young pulsars (Kargaltsev et al. 2006). Pulsa- 
tions were firmly detected only for three brighter pulsars, 
B1929+10, B0628-28, and B0950+08, with a hint of en- 
ergy dependence in the latter case, but their analysis was 
not very conclusive (e.g., Becker et al. 2004 and Zavlin 
& Pavlov 2004 came to quite different conclusions about 
the possible polar cap contribution to the emission from 
B0950+08). At this point, our working hypothesis is 
that the emission from old pulsars possibly includes both 
the magnetospheric and polar cap components, such that 
the magnetospheric component dominates at higher en- 
ergies, > 1-2 keV, and has about the same PL slope and 
efficiency as in young pulsars. However, we cannot rule 
out that old pulsars are purely magnetospheric emitters, 
with spectra softer, and efficiencies higher, than those of 
young pulsars. To check which of these two hypotheses is 

Newton observations of two other old pulsars, of which one, PSR 
B0834+06, was detected and showed properties similar to the pre- 
viously detected old pulsars. 
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Fig. 1. — Portion of the P-P diagram for pulsars from the 
ATNF Pulsar Catalogue (Manchester et al. 2005). The eight old, 
low-powered pulsars previously observed with Chandra and XMM- 
Newton are marked by crosses; J0108 is marked by the filled circle. 
The numbers 1 through 9 near the marked pulsars correspond to 
those in Fig. 9. Notice that J0108 is located close to the bottom- 
right boundary of the "swarm" of the pulsar population, generally 
known as the "pulsar death line". 

correct, not only deeper observations of X-ray brightest 
old pulsars are warranted, but also the number of ob- 
served pulsars should be increased. Particularly impor- 
tant would be X-ray observations of pulsars even older 
and less powerful than those in the currently observed 
sample. 

A very interesting target for such observations is PSR 
J0108-1431 (hereafter J0108), a very old, non-recycled 
pulsar (period P = 0.808 s, spin-down age r = P/2P = 
166 Myr, spin-down power E — 5.8 x 10 30 ergs s _1 , sur- 
face magnetic field B = 2.5 x 10 11 G), discovered by 
Tauris et al. (1994). Based on its dispersion measure, 
DM = 2.38 pc cur 3 (D'Amico et al. 1998), the smallest 
known for any radio pulsar, and the Galactic coordinates, 
/ = 140.9° and b = —76.8°, the distance to the pulsar is 
d = 130 pc, according to the model for Galactic electron 
density distribution by Taylor & Cordes (1993) (d = 180 
pc in the model by Cordes & Lazio 2002). J0108 is not 
only among the nearest NSs to Earth, but it also has the 
second lowest known radio luminosity, 0.15 mJy kpc 2 at 
400 MHz (for d= 130 pc), and its position in the P-P 
diagram is very close to the "pulsar death line" (see Fig. 
1). Mignani et al. (2003) carried out deep optical obser- 
vations of the J0108 field with the Very Large Telescope 
(VLT) and reported upper limits on its fluxes in several 
bands (e.g., V > 28.0, B > 28.6). 

Obviously, even a mere detection of such an extremely 
old, "almost dead" pulsar in X-rays would allow one to 
probe the pulsar properties in the regime where no other 
pulsar has been detected outside the radio range. There- 
fore, we observed J0108 with the Chandra X-ray Obser- 
vatory. We describe the observation and data analysis in 
§ 2, and discuss implications of our results in § 3. 

2. OBSERVATIONS AND DATA ANALYSIS 

The field of J0108 was observed by Chandra on 2007 
February 5 (MJD 54,136) using the Advanced CCD 
Imaging Spectrometer (ACIS). The observation was car- 



ried out in Very Faint mode, and the target was imaged 
on the ACIS-S3 chip, 7" off-axis. The detector was oper- 
ated in Full Frame mode, which provides time resolution 
of 3.24 s. The total scientific exposure time was 30,090 s. 
As no bad-time intervals (e.g., particle background flares) 
were found in the data, the whole exposure was used for 
the analysis. Data reduction was done with the Chandra 
Interactive Analysis Observations (CIAO) software (ver. 
3.4). For the analysis, we used the data reprocessed on 
2007 August 5 (ASCD ver. 7.6.11, CALDB ver. 3.4.1) 2 . 
We removed the pipeline pixel randomization and used 
the energy range of 0.3-8 keV to minimize the particle 
background. We used XSPEC ver. 11.3.2ag for spectral 
fits. 



Fig. 2. — Original ACIS image in the 0.3-8 keV band (pixels size 
0.492"). The 1.6" radius circle shows the source aperture from 
which the counts for spectral analysis were extracted. The cross 
shows the radio pulsar position from the ATNF Pulsar Catalogue 
(Manchester et al. 2005). 



2.1. Image 

Inspection of the ACIS image shows only one source in 
the vicinity of the radio pulsar position (see Fig. 2). To 
measure the coordinates of this pointlikc source, we used 
the CIAO tasks celldetect, wavdetect and vtdetect, which 
yielded the same source position within the centroiding 
uncertainty (w 0.05" for each of the coordinates): 

a = 01 h 08 m 08.354 s , S = -14°3l'50.38". (1) 

The astrometric uncertainty of this position can be esti- 
mated from the empirical distribution of the radial off- 
sets of X-ray positions (aspect solutions) from the ac- 
curately known celestial locations for a sample of point 
sources 3 . In particular, 68% of the 184 such sources 
imaged on the S3 chip within 2' from the optical axis 
have offsets smaller than 0.21" = -R0.687 and 90% of 
the sources have offsets < 0.44" = -Ro.90- For a two- 
dimensional Gaussian distribution with equal standard 
deviations a for the (uncorrelated) right ascension and 
declination offsets, we can estimate the standard devi- 
ation as u = R a {2 In [(1 — a) -1 ] } ! , where a is the 
probability to have a radial offset smaller than R a . This 

2 This reprocessing corrected for a systematic aspect offset of 
m 0.4", believed to be the result of a changed thermal environment 
in early 2007. 

3 See § 5.4 and Fig. 5.4 in the Chandra Proposers' Observatory 
Guide, ver. 10, at http://asc.harvard.edu/proposer/POG. 
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gives a = i?o. 68 /1.51 = 0.14" and a = R , 9 q/2A5 = 
0.205". The different values of a are likely due to the 
scarce statistics, and we can adopt the larger one as a 
conservative estimate. Adding the centroiding uncer- 
tainty in quadrature, we obtain the total uncertainty 
ir a — as — 0.21" for each of the coordinates. 

The absolute astrometric position of the target could 
be improved, in principle, by cross-correlation of the 
ACIS positions of field X-ray sources with the positions of 
their optical counterparts in catalogs with sufficiently ac- 
curate astrometry, such as the 2MASS catalog. However, 
we found only two X-ray sources with possible 2MASS 
counterparts. One of these sources, located 3.3' north- 
west of the target, has too few counts to measure its 
centroid with sufficient precision (centroiding errors are 
0.22" in right ascension and 0.11" in declination, while 
the corresponding Chandra — 2MASS offsets are 0.22" 
and —0.15"). As the other source, 3.7' southeast of the 
target, looks slightly extended (perhaps a double source) 
in X-rays, its centroid cannot be reliably measured (e.g., 
the positions measured with celldetect and wavdetect dif- 
fer by w 0.5"). Therefore, the field sources cannot be 
used for a reliable boresight correction. 

The measured X-ray position differs by 2.10" from the 
radio pulsar position in the ATNF Pulsar Catalogue 4 
(Manchester et al. 2005) derived from timing observa- 
tions of the pulsar (Hobbs et al. 2004) , and it differs by 
1.16" from the radio-interferometric position reported by 
Mignani et al. (2003). These differences significantly ex- 
ceed the absolute position uncertainty of Chandra source 
location and the reported uncertainties of the radio po- 
sitions (see Fig. 3). Therefore, the source could be either 
a field object (a star or an AGN) or the pulsar with an 
appreciable proper motion. Based on the lack of optical 
sources in the most conservative Chandra error circle, we 
will argue in § 2.2 that only the latter option is viable, 
and we will estimate the pulsar proper motion in § 3.1. 

2.2. Spectral Analysis 

Using the CIAO's psextract task, we extracted 53 
counts from the circular aperture of 1.6" radius (« 95% 
encircled energy radius) around the centroid position 
(Fig. 2). The background was taken from an annulus of 
9" < r < 18" around the pulsar. The background region 
has 57 counts, which corresponds to 0.60 counts when 
scaled to the source aperture. Thus, the source count 
rate is 1.74 ± 0.24 counts ks _1 , in the 1.6" aperture and 
0.3-8 keV energy range. In Figure 4 we show the ener- 
gies of all the 53 events detected, which are within the 
range 0.3 keV < E < 6keV, and the distribution of their 
arrival times over the duration of the observation. The 
distribution of the arrival times does not show statisti- 
cally significant deviation from the Poisson statistics, i.e. 
the source does not show flares, which could be seen in 
X-ray emission from a star or an AGN. 

The observed X-ray flux and its variance can be 
estimated as F = T^^2^ =1 / A), and a 2 F = 

T ^lJ2k=i( E k/ A k) 2 , where T cxp is the exposure time, 
N is the total number of detected events, Ek is the 
energy of fc-th event, and Ak — A(Ek) is the detec- 
tor effective area at energy Ek- For the source un- 

4 http:/ /www. atnf.csiro.au/research/pulsar/psrcat 



Fig. 3.— 11" X 10" VLT image of the J0108 field in the V 
band (1200 s exposure taken with the FORS1 instrument of the 
VLT-Antu in 2000 July 6-7; the reduced VLT image was kindly 
provided by R. Mignani). The position of the detected X-ray source 
is shown by the 0.21" radius Chandra error circle (68% confidence). 
The crosses 1.16" and 2.10" north-northwest of the X-ray source 
position correspond to the radio pulsar positions from Mignani et 
al. (2003) and Hobbs et al. (2004), respectively; the sizes of the 
cross arms are the uncertainties of radio positions. 
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Fig. 4. — Energies and arrival times for the 53 photons extracted 
from the 1.6" aperture around the pulsar. 



der consideration, the flux in the 0.3-8 keV band is 
^o.3-8kov = (9.2 ± 1.9) x 10~ 15 ergs cm" 2 s" 1 . This flux 
is much higher than the upper limits on the optical fluxes 
within the 5er Chandra error circle, Fy < 2.6 x 10~ 17 
and Fb < 2.4 x 10~ 17 ergs cm~ 2 s _1 in the V and B 
bands, respectively (Mignani et al. 2003). The large X- 
ray/optical flux ratio [e.g., Ry = -Fo.3-3.5 kev/Fv > 230, 
where F0.3_3.5kcV = (7-3 ± 1.2) x 10~ 15 ergs cm~ 2 s _1 ] 
means that the detected source cannot be an ordinary 
star or an AGN, the usual field sources at high Galactic 
latitudes, for which Ry < 0.3 and Ry < 50, respectively 
(see Fig. 1 in Maccacaro et al. 1988, and Table 1 and 
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TABLE 1 
Fits of the J0108 spectrum 
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a Hydrogcn column density in units of 10 2c> cm 2 . Fixed values 
are in square brackets. 

b Photon index for the PL fits or temperature for the BB fits, in 
eV. 

C PL normalization in units of 10" 6 photons cm" 2 s" 1 keV" 1 at 
1 keV or apparent emitting area in m 2 . 

d Unabsorbed flux in the 0.3—8 keV range, in units of 10 — 14 ergs 
cm" 2 s — 1 . 

°Best-fit C statistic value, for 526 bins used in the fit. 

-ray luminosity in the 0.3-8 keV band for the PL fits, or bolo- 
metric luminosity for the BB fits, in units of 10 2s ergs s — 1 . 

Fig. 1 in Stocke et al. 1991). The only known objects 
with R v > 100 are NSs (e.g., R v ~ 10 25 -10 4 for radio 
pulsars; see Zavlin & Pavlov 2004). Therefore, we con- 
clude that the source should be the X-ray counterpart 
of PSR J0108-1431, whose displacement from the previ- 
ously measured pulsar positions can be attributed to the 
proper motion of the pulsar. 

Because of the small number of counts, we used the C- 
statistic (Cash 1979) for spectral fitting, without binning 
the count spectrum. The results of spectral fits are pre- 
sented in Table 1. The errors of the fitting parameters 
in Table 1 and below in the text are quoted at the 68% 
confidence level for one interesting parameter, calculated 
with the XSPEC error command. 

First, we fit the spectrum with the absorbed PL model, 
with all the model parameters allowed to vary. Because 
of the small number of counts, the fit is poorly con- 
strained (see Fig. 5, top). Moreover, the best-fit hydro- 
gen column density obtained in this fit, iVn = 2.3 x 10 21 
cm" 2 , is unrealistically large, at least a factor of 10 larger 
than the total Galactic neutral hydrogen column density 
in the direction of J0108: N m w 1.8 x 10 20 cm" 2 (Dickey 
and Lockman 1990) or 2.1 x 10 20 cm" 2 (Kalberla et al. 
2005). The best- fit photon index obtained in this fit, 
T 3.4, is considerably higher than those usually ob- 
served from radio pulsars. 

To better constrain the model parameters, we have to 
fix the hydrogen column density. Using the pulsar's dis- 
persion measure (DM = 2.38 cm" 3 pc" 1 , which corre- 
sponds to the electron column density N e = 7.3 x 10 18 
cm" 2 ) and assuming the typical 10% ISM ionization, 
we obtain iV H = 7.3 x 10 19 cm" 2 . With the hydrogen 
column density fixed at this value, we obtain a smaller 
photon index, T « 2.2, and a lower unabsorbed flux, 

F o U 3-8kcV = L1 x 10_14 er S s cm ~ 2 s_1 ( see Fi S- 6 > to V, 
Fig. 7, and Table 1), corresponding to the isotropic lu- 
minosity L . 3 -8kcV ~ 2.1 x 10 28 d 2 30 ergs s" 1 (d 130 = 
<i/130pc). Changing the value of the fixed hydrogen 
column density within a reasonable range, Nn < Nm, 
changes the fitting parameters only slightly because the 
ISM absorption in the observed energy band is very low 
at such small Ah values (see the fit for A H = 2 x 10 20 
cm" 2 in Table 1). 
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Fig. 5. — 68% and 90% confidence contours in the Nn-T plane for 
the PL model (fop) and N^-kT plane for the BB model (bottom). 
The contours are obtained with the normalization parameters fitted 
at each point of the grids. 

The fit of the same count-rate spectrum with the ab- 
sorbed BB model with the hydrogen column density al- 
lowed to vary yields lower Ah values (Ah < 5 x 10 20 
cm" 2 , consistent with Ah = 7.3 x 10 19 cm" 2 estimated 
from the DM value; see Fig. 5, bottom, and Table 1). 
The fit with the column density fixed at Ah = 7.3 x 10 19 
cm" 2 gives approximately the same temperature kT w 
0.28 keV (T w 3.2 MK) and projected emitting area 
A± ~ 53d 2 30 m 2 as the fit with free Na (see Fig. 6, bot- 
tom, and Fig. 8) . This temperature and area correspond 
to the apparent radius R = (A±/n) 1 / 2 ~ 4.1di3o m and 
bolometric luminosity Lboi = iA±aT 4 ~ 1.3 x 10 28 d 2 30 
ergs s _1 . Similar to the PL fit, the fitting parameters 
very weakly depend on Nu (see Table 1) as long as it re- 
mains within reasonable limits. We also tried a PL+BB 
fit, but the fitting parameters could not be constrained 
because of the small number of counts. 

3. DISCUSSION 

3.1. Proper motion 

Based on the large X-ray/optical flux ratio, we have 
concluded in § 2.2 that the point source we detected with 
the Chandra ACIS is the X-ray counterpart of the J0108 
pulsar. As its X-ray position is different from those re- 
ported from earlier radio observations (see § 2.2), we can 
estimate the pulsar's proper motion. 
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Fig. 6. — Binned count-rate spectra for the data and the best- 
fit absorbed PL (top) and BB (bottom) models at fixed Nn = 
7.3 X 10 19 cm~ 2 (see Table 1). The counts in the data are binned 
to 5 counts per energy bin (except for the last bin that contains 3 
counts). The lower parts of the panels show the contributions of 
the spectral bins to the C-statistics (d — m is the difference between 
the data counts and model counts in an energy bin). 
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Fig. 7. — 68% and 90% confidence contours in the F-Af plane 
for the absorbed PL models with Nn values fixed at 2.3 X 10 21 
and 7.3 X 10 19 cm -2 . The PL normalization Af is in units of 
10 -6 photons keV -1 cm~ 2 s — 1 . The dashed lines are the lines of 
constant unabsorbed flux in the 0.3—8 keV band in units of 10~ 14 
ergs cm~ 2 s -1 . 

Among the previously measured pulsar coordinates, 
the most reliable ones, 

a = 01 h 08 m 08.317 s ± 0.010 s , S = -14°31'49.35" ± 0.35", 
were obtained in the 11.3 hr intcrferomctric observation 
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Fig. 8. — Confidence contours (68% and 90%) for the absorbed 
BB model with fixed N H = 7.3 X 10 19 cm" 2 . The BB normalization 
is expressed in terms of projected emitting area, in units of m 2 for 
d = 130 pc. The dashed lines are the lines of constant bolometric 
luminosity, in units of 10 28 ergs s -1 , assuming d = 130 pc. 

with the Australia Telescope Compact Array (ATCA) at 
epoch MJD 51,999 (Mignani et al. 2003).' Comparing 
these coordinates with those we measured, we obtain the 
proper motion 

fi a = 92±44mas/yr, fi s = -176 ± 70 mas/yr, (2) 



or 



H = 199 ± 65 mas/yr, P.A. = 152° ± 15° 



(3) 



for the total proper motion and the position angle 
(counted east of north). 

In addition to the interferometric position, timing posi- 
tions of J0108 have been reported by Tauris et al. (1994), 
DAmico et al. (1998), and Hobbs et al. (2004). The po- 
sition reported by Hobbs et al. (a = 01 h 08 m 08.30 s ± 
0.03 s , <5 = -14°31'48.4± 0.9"; mean epoch MJD 50,889 
[1998.3]), quoted in the current version of the ATNF Pul- 
sar Catalogue (Manchester et al. 2005), is supposed to be 
the most accurate because it used data acquired during 
the longest time span of 9.6 yr, versus 1 yr and 3 yr for 
the measurements by Tauris et al. and D'Amico et al., 
respectively. Moreover, the Hobbs et al. position sup- 
resedes the earlier ones because it includes the data used 
by Tauris et al. and D'Amico et al. 5 Using the Hobbs 
et al. position together with the the ATCA and Chandra 
positions, we find the proper motion 



fi a = 91 ± 40 mas/yr, fis 
/i = 210 ± 57 mas/yr, 



= -189 ±61 mas/yr, (4) 
P.A. = 154° ± 12°, (5) 



consistent with that given by equations (2) and (3). (We 
should note, however, that the weights of the ATCA- 
Chandra pair in this u, a and fis values, 0.76 and 0.73, 

5 We should note that the position reported by D'Amico ct al. 
is not consistent with the other two timing positions if its very 
small errors (0.005 s and 0.1" in R.A. and deck, respectively) are 
taken at face value. As these position errors are a factor of 10 
smaller than those of the Hobbs et al. position based on a factor 
of 3 longer observation time span, we suppose that the errors were 
strongly underestimated. Therefore, the upper limit on the proper 
motion, /Li < 82 mas/yr, obtained by Mignani ct al. (2003) from 
the comparison of the ATCA position with the position reported 
by D'Amico ct al., should be disregarded. 
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Fig. 9. — X-ray luminosities of nine old pulsars versus spin- 
down power. The numbers in parentheses correspond to the pulsars 
marked in Fig. 1. For PSR B0950+08 and PSR 1929+10, the 
diamond and the asterisk show the luminosities of the nonthermal 
and thermal components, respectively (as determined by Zavlin & 
Pavlov 2004 and Misanovic ct al. 2008, respectively). For the other 
pulsars, the luminosities were obtained from PL fits (see Fig. 5 in 
Kargaltsev et al. 2006). For the pulsars whose parallaxes have not 
been measured, we ascribe a factor of 2 uncertainty to the distances 
estimated from the pulsar's dispersion measure. 

respectively, are larger than those of the two other pairs 
because of the larger uncertainties of the timing posi- 
tion.) 

The inferred proper motion corresponds to the trans- 
verse pulsar velocity V± = (130±35)di3o km s -1 , some- 
what lower than the typical transverse velocity of 300- 
400 km s _1 . The proper motion in the Galactic coor- 
dinates is /ii w 140mas/yr, fj,b ~ — 155mas/yr, i.e., the 
pulsar is moving farther south from the Galactic plane. 

Interestingly, a recent re-analysis of the original VLT 
observations (epoch 2000.6) have shown a faint object 
(U = 26.4 ± 0.3, B « 27.9, V > 27.8), whose posi- 
tion (at the outskirts of the bright field galaxy) virtually 
coincides with that following from the backward extrap- 
olation of the proper motion to the epoch of the VLT 
observations (Mignani et al. 2008). Assuming this is the 
optical counterpart of J0108, its X-ray/optical flux ra- 
tios for the V and B bands (R v > 220, R B ~ 160) are 
consistent with typical values for radio pulsars (Zavlin & 
Pavlov 2004; Zharikov et al. 2006), while the relatively 
high flux in the U band might indicate thermal emission 
from the bulk of the NS surface with a brightness tem- 
perature ~9x 10 4 K, too low to be detected in X-rays 6 . 
The optical detection (if confirmed) provides additional 
support to the identification of the detected X-ray source 
with the JO 108 pulsar and to our estimate of the proper 
motion. 

3.2. Energetics and spectrum 

J0108 is the least powerful and the oldest among the or- 
dinary (nonrecycled) pulsars detected in the X-ray range. 
The previous record-holders have a factor of 15 higher 
spin-down power (B1133+16) and a factor of 9.5 smaller 

6 Such a relatively high surface temperature of the very old NS 
could be explained by (re)heating processes in the NS interiors (see 
Mignani ct al. 2008 for discussion and references). 



spin-down age (B0950+08). Therefore, our observation 
of J0108 has yielded the pulsar X-ray properties, such 
as the spectrum and luminosity, in the yet unexplored 
parameter domain, allowing us to examine the evolution 
of pulsar X-ray emission at very late stages of pulsar ac- 
tivity. 

The small number of source counts detected does not 
allow us to differentiate between the alternative spec- 
tral models for J0108. For the PL fit with the hydrogen 
column density fixed at TVr- = 7.3 x 10 19 cm~ 2 , we ob- 
tained the spectral slope F = 2.2 ± 0.2, somewhat softer 
than the typical values, T ~ 1-2, observed for young, 
more powerful pulsars (e.g., Gotthelf 2003). However, 
some old pulsars show even steeper spectra, T ~ 2-3 
(e.g. r = 2.5 ± 0.3 for PSR Bl 133+16; Kargaltsev et al. 
2006), suggesting that pulsar spectra might soften with 
increasing age (i.e., decreasing spin-down power). This 
softening can be ascribed to cither an evolution of the en- 
ergy spectrum of magnctosphcric particles or to a larger 
contribution of soft thermal emissison from polar caps in 
old, low-powered pulsars (e.g., Zavlin & Pavlov 2004). 

The luminosity found from the PL fit (e.g., 
£o.3-iokcV ~ 2.1 x 10 28 d 2 30 er S s s_1 ! see § 2-2 and Table 
1) is the lowest among the luminosities of ordinary pul- 
sars detected in X-rays. This is demonstrated in Figure 
9, where we plot the 1-10 keV luminosities of nine old, 
low-powered pulsars (E < 10 34 ergs s _1 ). Adding J0108 
to the sample makes it clear that the X-ray luminosity of 
old pulsars is correlated with the spin-down power in the 
same manner as has been observed for younger pulsars 
(i.e., Ex generally grows with E, albeit with a large scat- 
ter; e.g., Possenti et al. 2002; Li et al. 2008; Kargaltsev 
& Pavlov 2008). The X-ray efficiency, rj x = L x /E, can 
be estimated for J0108 as ?7o.3-8kcV ~ 4 x 10~ 3 d 2 30 . This 
efficiency is larger than the typical efficiencies, ~ 10 -5 - 
10~ 4 , of younger, more powerful pulsars (see e.g., Fig. 
5 in Kargaltsev & Pavlov 2008), but it is comparable 
with the efficiencies of old pulsars (see Fig. 9). This sup- 
ports the conjecture that, on average, old pulsars radiate 
a larger fraction of their spin-down power in the X-ray 
range than younger ones (Zharikov et al. 2006; Kargalt- 
sev et al. 2006). 

The higher X-ray efficiency could be attributed to the 
additional contribution of soft X-ray polar cap emission, 
in accordance with the softer X-ray spectra of old pul- 
sars. This hypothesis could be tested by deep X-ray ob- 
servations that would allow phase-resolved spectroscopy 
and energy-resolved timing (because the polar cap radi- 
ation, which dominates at lower energies, should show 
shallower pulsations than the magnctosphcric radiation; 
e.g., Zavlin & Pavlov 2004). Alternatively, one could 
speculate that the broadband magnetospheric spectrum 
shifts toward lower energies with decreasing E, so that 
the soft X-ray efficiency increases at the expense of de- 
creasing hard X-ray or soft 7-ray efficiency. Observations 
at higher energies with future missions could verify this 
hypothesis. 

As we have mentioned in §2.2, the observed X-ray spec- 
trum of J0108 can be fitted with a BB model, which 
might be interpreted as thermal radiation from heated 
polar caps. The best-fit BB temperature, 3.2 MK, ob- 
tained from this fit, is within the 1.7-3.5 MK range found 
from BB (or PL+BB) fits for other old pulsars (Zavlin 
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& Pavlov 2004; Kargaltsev et al. 2006; Misanovic et al. 
2008; Gil et al. 2008). The projected emitting area, 
A± ~ 50d 2 30 m 2 , is, however, exremely small. Not 
only it is considerably smaller than the smallest area 
reported (~ 240 m 2 for PSR B0834+06), but also the 
ratio, ~ 6 x 10 — 4 df 30 , of this area to the conventional 
polar cap area, A pc = 2-k 2 R^ s /(cP) ~ 8.7 x 10 4 m 2 , 
is smaller than for any other pulsar. We should bear 
in mind, however, that both the temperature and the 
area are not only highly uncertain because of the strong 
correlation of these parameters (see Fig. 8), but they are 
also strongly dependent on the model of thermal emission 
(e.g., fits with hydrogen NS atmosphere models usually 
give a factor of 2 lower effective temperature and a fac- 
tor of 10-100 larger area; Pavlov et al. 1995; Zavlin & 
Pavlov 2004). The large difference between the observed 
and conventional areas can be also explained by the as- 
sumption that the non-dipolar component of the mag- 
netic field at the NS surface is much stronger than the 
dipolar component inferred from the pulsar spindown, in 
accordance with the partially screened gap model for the 
inner acceleration region above the pulsar polar cap (Gil 
et al. 2003). Finally, the relationship between the actual 
emitting area and A± depends on the orientation of the 
magnetic and spin axes, and on the General Relativity 
effects (see discussion in Pavlov et al. 2007). 

The apparent bolometric luminosity, Lboi ~ 1-3 x 
10 28 <i 2 3g ergs s" 1 , which is less mo del- dependent than 
Aj_ and T (Zavlin & Pavlov 2004), is also smaller than 
those of the other old pulsars; however, the correspond- 
ing "polar cap efficiency", n pc = L^ \/E ~2x 10~ 3 d 2 30 , 
is comparable to (or even larger than) those reported 
for other pulsars. Such an efficiency exceeds the pre- 
dictions of the polar cap heating models by Harding & 
Muslimov (2001, 2002) 7 , but it is in a reasonable agree- 
ment with the partially screened gap model by Gil et al. 
(2003, 2008) . We should not forget, however, that this 
r/ pc is only an upper limit because the observed spectrum 
can contain a nonthermal component. 

7 Since the curvature radiation cannot induce a pair cascade 
in the low-powered J0108, only the less efficient inverse Compton 
scattering (ICS) cascade can be responsible for polar cap heating 
in this pulsar. The Lbol and ?7pc estimated from our observation 
of J0108 are close to the analytical estimates for upper limits on 



3.3. Summary 

Thanks to the high sensitivity of the Chandra ACIS, 
we have detected J0108, the oldest and the least pow- 
erful pulsar ever observed in X-rays. This observation 
has shown that even very old pulsars can be rather ef- 
ficient X-ray emitters, even more efficient than young 
ones. It seems likely that the observed emission contains 
both magnetospheric and thermal components. To sep- 
arate these components and understand the nature of 
the X-ray emission, a deep XMM-Newton observation, 
which would provide not only high sensitivity but also 
a sufficient time resolution, would be particularly use- 
ful. To evaluate the X-ray luminosity and efficiency of 
JO 108 more accurately, the distance to the pulsar should 
be determined from annual parallax measurements. 

Comparing the Chandra position of JO 108 with the 
previously measured radio positions, we have been able 
to estimate the pulsar's proper motion of about 200 
mas/yr. To determine the proper motion with a higher 
accuracy, radio observations with ATCA and/or VLA are 
required. Alternatively, the proper motion can be accu- 
rately measured with new Chandra observations in a few 
years from now. 

The detection of the candidate optical counterpart of 
JO 108 (Mignani et al. 2008) was hampered by the prox- 
imity of its sky position to the bright field galaxy (see Fig. 
3). Our estimate of the proper motion suggests that it 
has moved away from the galaxy and could be firmly de- 
tected in the near future. Measuring its optical-UV spec- 
trum would be particularly important for understanding 
the thermal evolution and magnetospheric emission of 
old NSs. 



We thank Roberto Mignani for providing the reduced 
VLT images of the pulsar field and useful discussions. 
This work was partially supported by Chandra award 
SV4-74018. 

these quantities for the resonant ICS heating (eqs. [60] and [64] 
in Harding & Mulsimov 2002). However, an extrapolation of the 
numerical calculations of these authors into the J0108 parameter 
domain seems to predict lower values. 
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